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Abstract:The exploration of spin symmetry (SS) in nuclear physics has been instrumental in identifying 
atomic nucleus structures. In this study, we solve the Dirac equation from the relativistic mean field (RMF) 
in complex momentum representation. We investigated SS and its breaking in single-particle resonant states 
within deformed nuclei, with a focus on the illustrative nucleus ‘°*Er. This was the initial discovery of a reso- 
nant spin doublet in a deformed nucleus, with the expectation of the SS approaching the continuum threshold. 
With increasing single-particle energy, the splitting of the resonant spin doublets widened significantly. This es- 
calating splitting implies diminishing adherence to the SS, indicating a departure from the expected behavior as 
the energy levels increase. We also analyzed the width of the resonant states, showing that lower orbital angular 
momentum resonances possess shorter decay times and that SS is preserved within broad resonant doublets, as 
opposed to narrow resonant doublets. Comparing the radial density of the upper components for the bound-state 
and resonant-state doublets, it becomes evident that while SS is well-preserved in the bound states, it deterio- 
rates in the resonant states. The impact of nuclear deformation (G2) on SS was examined, demonstrating that an 
increase in 32 resulted in higher energy and width splitting in the resonant spin doublets, which is attributed to 
increased component mixing. Furthermore, the sensitivity of spin doublets to various potential parameters such 
as surface diffuseness (a), radius (R), and depth (Xo) is discussed, emphasizing the role of these parameters 
in SS. This study provides valuable insights into the behavior of spin doublets in deformed nuclei and their 


interplay with the nuclear structure, thereby advancing our understanding of SS in the resonance state. 
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I. INTRODUCTION 


Substantial spin symmetry (SS) breaking between spin 
doublets (n, l, 7 = | + 1/2) is one of the most important 
concepts in nuclear structure analysis, and has been exten- 
sively discussed in scientific literature [1-4]. SS is charac- 
terized by a quasi-simple single-nucleon doublet defined by 
quantum numbers (n, l, 7 = l + 1/2). SS not only facili- 
tates the understanding of the nuclear phantom number of the 
atom, but also forms a cornerstone of the structure of the nu- 
clear shell, providing an important framework for interpret- 
ing the spatial arrangement of nuclei [5, 6]. In the context 
of atomic nuclei, SS and spin-orbit (SO) splitting phenomena 
are crucial for understanding the structure of atomic nuclei. 
The interaction between the intrinsic spin and orbital angu- 
lar momentum leads to distinct energy levels for nucleons in 
the nucleus, contributing to the rich and complex behavior of 
nuclear systems [7, 8]. 

The discovery of SS has inspired extensive exploration of 
its origins, enhancing the understanding of nuclear structures. 
The extension of the shell model by Nilsson et al. [9, 10] 
provides a crucial framework not only for characterizing de- 
formed nuclei, but also for understanding phenomena related 
to nuclear rotation. Haxel et al. established the crucial sig- 
nificance of the SO potential [8], and Mayer et al. [6] demon- 
strated notable state splitting associated with an elevated or- 
bital angular momentum. When integrated with the mean- 
field potential energy, such as the Woods-Saxon potential or 
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harmonic oscillator potential, the classical magic number in 
nuclear physics can be efficiently reproduced using the forced 
SO potential. This amalgamation of SS, extended-shell mod- 
els, and SO potential collectively establishes a fundamental 
comprehension of the nuclear structure, augmenting our pro- 
ficiency in forecasting and elucidating various nuclear phe- 
nomena. 


Smith et al.[11] and Bell et al.[12] observed that the antin- 
ucleon spectrum approximates that of SS. Serot and Walecka 
successfully predicted SO splitting using relativistic mean 
field (RMF) theory [13]. Ginocchio [14] demonstrated SS 
and U(3) symmetry in the Dirac Hamiltonian under scalar 
and vector harmonic oscillator potentials (V(r) — S(r) = 0) 
within the RMF framework. Meng et al. , [15] demonstrated 
the exactness of SS when the condition d(V (r) — S(r))/dr = 
0 is satisfied. In the process of developing the RMF the- 
ory based on the Dirac Woods-Saxon field theory proposed 
by Zhou et al. [16, 17], they explored SS and pseudospin 
in nucleon and antinucleon spectra within RMF theory. Al- 
haidari et al.[18] investigated the physical interpretation of 
the three-dimensional Dirac equation under SS conditions 
[V(r) = S(r)]. Application of similarity renormalization 
group theory, as demonstrated by Guo et al. [19], pro- 
vided insight into the SS of the Dirac Hamiltonian under ax- 
ially deformed scalar and vector potentials. In addition, ref- 
erences [2] and related studies provide additional information 
on relativistic symmetry. This combination of experimental 
observations and theoretical frameworks contribute compre- 
hensively to the evolving understanding of SS in antinucleon 
spectra in the field of nuclear physics. 


In recent years, the exploration of single-particle reso- 
nant states, particularly in exotic nuclei with unconventional 
neutron-to-proton (N/Z) ratios, has gained significant atten- 
tion because of the discovery of novel phenomena [20-29]. 


Exotic systems, which often feature weakly bound or un- 
bound nuclei, are open quantum many-body systems that play 
key roles in the continuum [30, 31]. In these nuclei, the Fermi 
surface is positioned near the continuum threshold, rendering 
the valence nucleons prone to scattering into the continuum, 
which is influenced by the pairing correlations. Hence, gain- 
ing a comprehensive understanding of the contribution of a 
continuum is of utmost importance [32-47]. Therefore, in- 
vestigating the SS within single-particle resonant states be- 
comes significant. We reviewed key studies that employed 
diverse theoretical approaches. Zhang et al. [48] confirmed 
the similarity between the components of the Dirac wave- 
function of resonant-spin double-dipole state. Xu et al.[50] 
explored solutions to the Dirac equation with various poten- 
tials, establishing correlations between SO splitting and po- 
tential parameters, building on the work in Ref.[49]. Li et 
al. [51] applied the coupling channel method to the Dirac 
equation with a quadrupole-deformation Wood-Saxon poten- 
tial. Shi et al.[52], provided insights into resonances in the 
Yukawa potential using the complex scaling method, reveal- 
ing connections between SS quality and potential parameters. 
Xu et al.[53] incorporated SO coupling, resonance contribu- 
tion, and pairing correlation to describe deformed nuclei. Re- 
cent endeavors include those of Sun et al. [54], where the 
Green’s function method was used to investigate the SS in 
single-particle resonance states, and Shi et al. [55] ex- 
plored SS in Pb isotopes through the RMF-CMR theory to 
study spin properties in nuclear systems. These studies are 
aimed at identifying the correlations between the SS quality 
and a spectrum of parameters; the chronological progress un- 
derscores the continuous and evolving exploration of nuclear 
structures. 


As emphasized in Ref.[56], understanding SS in resonant 
states within deformed systems remains an unsolved prob- 
lem. This study employs a complex momentum representa- 
tion (CMR) approach, which offers several advantages over 
other methods. For example, CMR allows the equitable treat- 
ment of both bound and resonant states, enabling the simulta- 
neous determination of narrow and broad resonances. Recent 
developments in the CMR approach are discussed in detail in 
Refs. [57—60]. This study focused on scrutinizing the SS and 
its breaking into realistically deformed nuclei using the CMR 
method. The analysis considered the energies and widths 
of the spin doublets and the similarity of the upper compo- 
nents of the Dirac spin measurements. Furthermore, corre- 
lations between SO splitting in resonant-state doublets and 
the parameters of the Wood—Saxon potential have been ex- 
plored [50, 61, 62]. Section 2 outlines the theoretical frame- 
work, and Section 3 provides numerical details and results, 
culminating in a comprehensive summary in Section 4. This 
study aims to enhance the understanding of SS in the res- 
onant states of deformed nuclei, with potential implications 
for wider applications in nuclear physics. 


Il. FORMALISM 


From the RMF theory [63, 64], the Dirac equation is ob- 
tained as 


H |p) = [a- p+ (M +S) +V] W) =el), 0 


where @ and 8 are the Dirac matrices, W and p are the nu- 
cleon mass and momentum, respectively, and S and V are the 
attractive scalar and repulsive vector potentials, respectively. 
E = e — M represents the single particle energy, and W is the 
wavefunction. 

To obtain the resonant states, the Dirac equation (1) was 
expressed as follows: 


pe (zla i) ab (K) Sg (K) , (2) 


For axially deformed nuclei, the third component m; of the 
total angular momentum j and parity 7 are suitable quantum 
numbers. The Dirac spinor can be expanded as 


=m ()-C(p woes) © 


where the angular partial wave function is a two-dimensional 
spinor $1jm,(Qk) = Ym, (lmZms|IM;)Vim(ax)Xm, and 
(=2j-1. 

Substituting the Dirac wave function given in Eq. (3) into 
Eq.(2), the radial Dirac equation is transformed, resulting in 
a set of coupled-channel equations. These coupled-channel 
equations are then formulated in matrix form, which is sub- 
sequently diagonalized to obtain the solutions. This method- 
ology allows for the simultaneous determination of all bound 
and resonant states. Additional details can be found in [57]. 


Til. NUMERICAL DETAILS AND RESULTS 


This study aims to understand the complex behavior of 
nuclear particles by investigating spin doublets in resonant 
states. It seeks to unveil insights into the structure and prop- 
erties of deformed nuclei, with a focus on the energy levels, 
wave functions, and spin symmetry (SS) in resonant states. 
The nuclear potential is defined as follows: 


U(r) = Uof (r) — B2Uok (r) Yao (0, p), 


where 2 is the quadrupole deformation parameter. In our 
study, we adopt a potential, characterized by the functions 


f(r) = —ke and k(r) = aftr) Here, U(r) denotes both 
lte a 


the vector and scalar potentials, a denotes the diffusivity, and 
R is the radius. The adoption of a potential resembling the 
Woods-Saxon shape underscores the relevance of such math- 
ematical representations for capturing the intricacies of nu- 
cleon interactions within the nuclear environment, making 
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Fig. 1. The single neutron energies are determined for the spin dou- 
blets Q = A + 1/2[N,nz, A] in '®8Er with 82 = 0.34. Here, 
N represents the main quantum number of the harmonic oscillator, 
n3 denotes the number of quanta for oscillations along the symme- 
try axis (in the z-direction) and A and Q respectively represent the 
components of orbital and total angular momentum projected along 
the symmetry axis. 


them realistic enough for application to nuclei. The param- 
eters of this potential can be determined by fitting it to the 
self-consistent microscopic potential obtained from the RMF 
calculations using the NL3 effective interaction, where the 
mean field potential is © = V +S and A = V — S. For !®Er, 
the specific parameter values are a = 0.66 fm, R = 6.74 fm, 
No = —62.8 MeV, and Avg = 637.4 MeV. Our exploration of 
the properties of spin doublets in deformed nuclei will enable 
us to discuss and analyze SS in a deformed context, thereby 
providing valuable insights into its nuclear structure. 

Using these parameters, we solved the Dirac equation us- 
ing the CMR method to obtain the single-particle energy 
levels of both the bound and resonant states in ‘°°Er with 
Bo = 0.34. In Fig. 1, the spin doublets for the resonant states 
are presented along with the bound states. The panels, or- 
dered from left to right, illustrate varying A values ranging 
from 1 to 4. Remarkably, our analysis revealed 13 pairs of 
spin doublets for single-neutron resonant states. This obser- 
vation marks a significant milestone as it is the first instance 
of identifying resonant spin doublets in deformed nuclei. Fig- 
ure | illustrates that a good SS is anticipated in the vicinity of 
continuous thresholds, because in this region, the condition 
for the exact SS, d(V(r) — S(r))/dr = 0, is more likely to be 
fulfilled. Certain spin doublets, such as [871]1/2, [871]3/2, 
[743]5/2, [743]7/2, [734]7/2, and [734]9/2, present desirable 
SS near the energy threshold. However, as the single-particle 
energy increases, the splitting of the doublets increased, and 
the quality of the SS deteriorated when the spin doublets 
moved further from the continuous threshold. For all bound 
and resonant spin heavy states, it is noteworthy that the en- 
ergy state of the spin-down is consistently higher than that of 
the spin-up. 

The resonant states were characterized by their energies 


TABLE 1. The four resonant spin doublets in 165Er exhibit energies 
and widths measured in MeV with G2 = 0.34. 


A [N, n3, AJQ E,[MeV] T[MeV] 
1 [871]1/2 2.111 24.633 
[871]3/2 1.971 23.359 
2 [862]3/2 2.318 24.700 
[862]5/2 2.052 23.700 
3 [613]5/2 4.604 2.183 
[613]7/2 3.301 0.776 
4 [622]7/2 3.090 3.512 
[622]9/2 0.790 0.165 


and widths. Nuclei situated far from the stability valley ex- 
hibit a Fermi surface closely aligned with the continuum, ren- 
dering the valence nucleons prone to scattering. The lifetimes 
of the resonant states are particularly significant when assess- 
ing the potential formation of exotic phenomena. A strong 
correlation exists between the width of the resonance state 
and its duration. Notably, single-particle states with quantum 
numbers (n, l, 7 = l — 1/2) are generally wider than their 
spin-doublet counterparts (n, l, 7 = 1 + 1/2). These findings 
imply that resonant states with a lower orbital angular mo- 
mentum have shorter decay times. The lower centrifugal bar- 
rier associated with the lower orbital angular momentum con- 
tributed to the larger resonant widths. The neutron spin dou- 
blets near the continuum threshold are presented in Table 1. 
The energy and width of the wide resonance spin doublets are 
very similar, for example, [871]1/2,[871]3/2, [862]3/2, and 
[862]5/2. For the narrow resonant states, differences emerged 
in the energies of [613]5/2,[613]7/2, [622]3/2, and [622]5/2, 
with relatively significant discrepancies noted in the widths of 
these two pairs. These observations indicate the preservation 
of SS in broad resonant doublets. 


To underscore the significance of SS assessment, as em- 
phasized in Ref. [67], it is crucial to evaluate SS by com- 
paring the upper components of Dirac spinors. In a previous 
study [55] that focused on spherical nuclei, the resonant states 
were observed to share similarities with their upper compo- 
nents. However, challenges arise in deformed nuclei where 
spin doublet states exhibit wave functions with distinct spher- 
ical components, making direct comparisons between differ- 
ent doublets difficult. The radial density distribution compris- 
ing diverse spherical components is significantly correlated 
with the wave function. The Dirac spinor comprises an upper 
Dirac spinor f(r) and a lower Dirac spinor g(r). By calcu- 
lating the radial density of the upper component for the spin 
doublets, we can utilize pm; (r) = Ly [fi (r) f" (r)], 
which provides a means of assessing similarities. In Fig. 2, 
subfigure (a) clearly demonstrates the similarity between the 
radial densities of the upper components of the bound dou- 
blets [512]3/2 and [512]5/2. Moreover, subfigure (b) presents 
the real part of the radial density for the resonant doublets 
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Fig. 2. The radial density distributions of the spin doublets in 1°*Er 
with 62 = 0.34. The upper components of the radial density distri- 
butions for the bound spin doublet ({512]3/2, [512]5/2) and the real 
part of upper components for the resonant spin doublet ([613]5/2, 
[615]7/2) are shown in the subfigures (a) and (b). 


[613]5/2 and [613]7/2. A comparison with the bound dou- 
blets reveals similarities and highlights unique features in the 
radial density of the resonant spin doublets. The bound states 
exhibit well-maintained SS, whereas the quality of SS in the 
resonant doublets is somewhat diminished. The radial den- 
sity variations of the upper components in relation to their 
amplitudes within the spin doublets appear to be due to the 
effect of deformation and multiple spherical components in 
the deformed nuclei wavefunction. 

To gain a deeper understanding of the SS quality of the de- 
formed nuclear resonance states, we investigated the impact 
of the deformation parameter {2 on the SS. The results are 
presented in Fig. 3 show case plots illustrating the energy 
splitting and width splitting between the spin doublets. Sev- 
eral resonant and bound spin doublets exhibited variations in 
Bo. Fig. 3 clearly distinguishes the bound (solid lines) and 
resonant (dashed lines) spin doublets. The spherical labels 
at the 62 = 0 position provide a reference for comparison. 
In the range 62 = 0 to fg = 0.60, the energy splitting of 
the bound spin doublets exhibited fluctuations but remained 
predominantly unaffected by the quadrupole deformation pa- 
rameter $2, which is consistent with the findings in Refs. 
[19, 68]. For resonant spin doublets, the deformation pa- 
rameter (2 significantly influences the energy splitting. As 
{2 increases, energy splitting within the most resonant spin 
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Fig. 3. The calculated single-neutron energy and width splittings as 
a function of the quadrupole deformation 62, where the bound and 
resonant doublets are marked respectively by the solid and dashed 
lines. 


doublet also increases. The width splitting of the resonant 
spin doublets ([622]3/2, [622]5/2) and ([613]5/2, [613]7/2) 
increased with the deformation parameter (2, whereas the 
width splitting of the other two pairs of doublets initially de- 
creased and then increased as 8 increased. Thus, from AE 
and AT, it can be seen that the energy and width splitting be- 
tween the spin-partner states of the resonance state increases, 
and thus, the quality of the spin symmetry deteriorates as the 
deformation parameter increases. 

To examine the deformation parameter 6 and its impact 
on the increased spin splitting of the resonant spin doublets 
in 168ẸEr, we investigated the occupation probabilities of the 
major components within the doublet states [613]5/2 and 
[613]7/2, as presented in Fig. 4. Substantial changes in 
the occupation probabilities of the various components were 
observed across a range of deformations. At the [613]5/2 
level, the components ds/2, 97/2, 99/2, and gi1/2 display 
notable occupation probabilities. At the [613]7/2 level, the 
components 99/2, 97/2) 11/2, and 713/2 become prominent. 
The contributions of the other components were negligible. 
At 2 = 0, where the spherical configuration prevails, the 
[613]5/2 level corresponds to the degenerate spherical con- 
figuration of g7/2, whereas the [613]7/2 level is character- 
ized by the degenerate spherical configurations of gg/2 and 
97/2. These spherical configurations, specifically 97/2 and 
99/2, represent exemplary spin doublets characterized by a 
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Fig. 4. Occupation probabilities of major configurations as a func- 
tion of G2 for the resonant spin doublets [613]5/2 and [613]7/2 in 
168 

Er. 


minimal energy splitting of only 1.45 MeV and width split- 
ting of just 0.66 MeV. However, as the deformation parameter 
{2 increases, we observe a decrease in the occupation prob- 
abilities of g7/2 and gg/2 accompanied by an increase in the 
probabilities of the other components. At G2 = 0.60, the oc- 
cupation probabilities of g7/2 and gg/2 were reduced to 57% 
and 71%, respectively, resulting in an increase in the spin en- 
ergy splitting to 1.67 MeV and width splitting to 2.36 MeV. 
Therefore, it can be inferred that as the deformation increases, 
the mixing of components leads to deterioration of the SS. 


To elucidate the dynamics of the evolution of SS within the 
resonance state, we studied the correlation between spin split- 
ting and the intrinsic parameters of the Woods-Saxon poten- 
tial. Fig. 5 provides a visual representation of the relationship, 
highlighting the interplay between the between spin splitting 
and various parameters associated with the Woods-Saxon po- 
tential. The graphical representation distinguishes between 
resonant and bound doublets using open symbols for resonant 
doublets and filled symbols for bound doublets. We investi- 
gated the effects of surface diffuseness on energy and width 
splitting parameter while maintaining Xo and R constant. The 
visualization in Fig. 5 subfigure (a) illustrates the effect of the 
surface diffusivity parameter a. As a increases from 0.47 fm, 


the energy splitting of most spin doublets decreases, indicat- 
ing an improved quality of SS. In addition, there is a con- 
sistent decrease in the width splitting of all the spin doublets 
with increasing a. Notably, the sign of the spin splitting for 
the spin doublets does not reverse with this parameter change. 
Analyzing single-particle states with quantum numbers (n, l, 
j = l — 1/2), it was observed that their widths were consis- 
tently larger than their spin doublets (n, l, j = l + 1/2), sug- 
gesting longer decay times for the resonant states with higher 
orbital angular momentum. This phenomenon is attributed to 
the centrifugal barrier effect, in which a lower orbital angular 
momentum results in a lower centrifugal barrier, leading to 
larger resonance widths. In Fig. 5 subfigure (b), we present 
information on the effect of the radius (R) on energy and 
width splitting. Variations in these splittings are plotted as 
a function of R while holding the other parameters constant. 
Observations with increasing R (starting at R=6.54 fm) show 
that the energy splitting of bound spin doublets decreases, 
whereas that of the resonance states increases. In terms of 
the width-splitting trends, most spin doublets experienced a 
decrease as R increased, with the exception of [732]3/2 and 
[732]5/2. Particularly, there is a consistent width discrepancy 
in the single-particle states, where the widths of all quantum 
number states (n, l, 7 = l — 1/2) are consistently larger than 
their spin doublets (n, l, j = l + 1/2). This implies that 
even if the energies of the spin doublets are perfectly simple 
and fused, their decay times are different. The importance 
of these results are underlined by the isotope comparison, as 
the dependence on R is particularly striking when contrasting 
different isotopes, owing to noticeable changes in the radius 
of the mean-field potential. Finally, we investigated the effect 
of varying Xo while keeping the values of a and R constant, 
as shown in subfigure (c). Observations with varying |Xo 
show that as |X| decreases, the energy splitting of the bound 
spin doublets remains almost unchanged. The energy split- 
ting of the resonant spin doublet states decreased with de- 
creasing |Xo|. The splitting width generally decreased with 
increasing |Xio|. The width of the single-particle resonance 
states, which is an indication of the decay time, shows that a 
higher potential depth increases the stability of the resonance 
states. These results highlight the sensitivity of SS to parame- 
ter variations within the Woods-Saxon potential, particularly 
in terms of the energy and width splittings for the bound and 
resonant spin doublets. 


IV. SUMMARY 


This investigation focuses on the spin doublets within the 
single-particle resonant states of an axially deformed nucleus 
by employing the Dirac equation within the RMF framework. 
The nuclear potential is defined by the depth (Uo), deforma- 
tion parameter (82), and spherical harmonic functions. Us- 
ing !°°Er as an illustrative example, we calculated the self- 
consistent microscopic potentials and single-neutron energy 
levels for both the bound and resonant states. Thirteen pairs 
of spin doublets were identified in 188Er, exhibiting varying 
SS. These resonant spin doublets displayed varying degrees 
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Spin energy and width splittings as a function of every potential parameter for the resonant states shown in Fig. 1. Here, the 


quadrupole deformation 82 = 0.34. The data corresponding to the variables a, R, and Xo are respectively displayed in (a), (b), and (c) with 
the other parameters fixed to R = 6.74 fm, No = —62.8 MeV in (a), a = 0.67 fm, Xo = —62.8 MeV (b), and a = 0.67 fm, R = 6.74 fm in (c). 


of SS, with some approaching a good SS near continuous 
thresholds, whereas the quality of the SS deteriorated as they 
moved away from the threshold. A comparison of the ener- 
gies and widths of these resonant spin dual states indicated 
that the spin-up state possessed higher energy than that of the 
spin-down state. 

The radial density distributions of the upper components of 
the resonant states were examined, which indicated similari- 
ties and distinctive features between the different doublets. 
SS is well preserved in bound states but deteriorates in res- 
onant states owing to spin breaking. We also conducted a 
specific investigation of the impact of the quadrupole defor- 
mation parameter (82) on the SS. As 62 increased, there was a 
notable increase in the energy and width splitting of the most 
resonant spin doublets. In addition, the SS of the bound spin 
doublets remained relatively unaffected. Furthermore, an in- 
depth analysis of the occupation probabilities within the res- 
onant doublets was conducted. Observable changes in these 
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probabilities were observed as the deformation parameter {2 
increased, contributing significantly to the deterioration of SS 
in the resonant states. 


Probing spin splitting and its relationship by varying the 
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